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The cross sections of the ^He(a,a)''He and ^He(Q,7)^Be reactions are studied at low energies using 
a simple two-body model in combination with a double-folding potential. At very low energies the 
capture cross section is dominated by direct s-wave capture. Ifowever, at energies of several MeV 
the d-wave contribution increases, and the theoretical capture cross section depends sensitively on 
the strength of the L — 2 potential. Whereas the description of the L = 2 elastic phase shift requires 
a relatively weak potential strength, recently measured capture data can only be described with a 
significantly enhanced L = 2 potential. A simultaneous description of the new experimental capture 
data and the elastic phase shifts is not possible within this model. Because of the dominating 
extranuclear capture, this conclusion holds in general for most theoretical models. 

PACS numbers: 25.55.-e,25.55.Ci,24.50.-|-g,25.40.Lw 



I. INTRODUCTION 

The ^He(a,7)^Be capture reaction is a key reaction 
in nuclear astrophysics. The flux of solar neutrinos 
at higher energies depends on the branching between 
the '^He('^He,2p)a and '^He(a,7)'^Be reactions, and in 
big-bang nucleosynthesis ''Li can be produced via the 
'^He(a,7)^Be capture reaction and subsequent electron 
capture of ^Be While the energy range of big-bang 
nucleosynthesis is covered by experimental data, the ex- 
periments are approaching the Gamow window around 
Eq — 22keV of the ^He(Q!,7)^Be reaction in the sun, but 
still have not reached Eq. 

Very recently, it has been attempted successfully to 
extend the measured energy range of the 3He(a,7)^Be 
reaction to higher energies up to about 3.2 MeV 0. 
From these new experimental data the energy depen- 
dence of the '^He(a,7)^Be capture reaction can be ex- 
tracted and compared to theoretical predictions. A bet- 
ter theoretical understanding of the energy dependence 
will help to reduce the uncertainties of the extrapola- 
tion of the cross section down to the Gamow window 
at the temperature of the interior of our sun. Several 
further experimental data sets exist at lower energies 
I, i, i, [6, 7, 8, 9, 10, [m, [13, [H, Ei, [H, p that 
are summarized in two recent compilations 18|, Il9l | . 

Besides theoretical calculations shown in the experi- 
mental papers I, i, s i 0, a a 0, [n m, [Hi 0, 

[Tsl . [T^ . \n\ , various theoretical studies have been devoted 
to the analysis of the '^He(a,7)^Be ca ptu re cross section 

fMMMM M, M [23, In, [H M [m, [H M, M 
|36| . ISTl |38| . 139( 1 . Additionally, a review on the status 
of the '^He(Q!,7)^Be reaction is given in [ioj . 

Theoretical models can be significantly constrained by 
the request that elastic scattering data have to be de- 
scribed simultaneously with the capture data. Phase 
shifts have been derived from elastic scattering angular 



distributions in [4l|, [42, [43|, and angular distributions 
or excitation functions at low energies are reported in 

[3ii[il[ii, [li- 
lt will be shown that the capture cross section of the 
■^He(Q;,7)^Be reaction at low energies is dominated by El 
transitions from incoming s- and c?-waves to bound p- 
states in ^Be. Consequently, a precise description of the 
s-wave and d-wave phase shifts is a prerequisite for the 
calculation of the ^He(Q:,7)^Be capture cross section. 

This article is organized as follows. In Sec. [11] a brief 
description of the direct capture model is provided. Gen- 
eral remarks on the applicability and reliability of this 
model and other theoretical calculations are given in 
Sec. mil Results for the ■^IIe(a,a)^He elastic scattering 
and ^He(a,7)^Be capture cross sections are presented in 
Sec. IIVI and theoretical uncertainties are carefully ana- 
lyzed. Conclusions are drawn in Sec. |Vl and finally a 
brief summary is given in Sec. IVII 

All energies E are given in the center-of-mass system 
throughout this paper except explicitly noted. Excitation 
energies E* in "^Be are related to E hy E* — E + Q with 
the Q-value of the •^IIe(Q:,7)''Be reaction Q — 1586.6 keV 
lil. 



II. THE DIRECT CAPTURE MODEL 

The direct capture (DC) model is a simple but pow- 
erful model to calculate capture cross sections between 
light nuclei. The full formalism is given explicitly in [isj . 
and its application together with systematic folding po- 
tentials is described in detail in our previous work on the 
•^IIe(Q!,a)^He and ^He(a,7)^Be reactions [3l|. Here I only 
repeat some essential features of the DC model that are 
important for the following discussion. 

The theoretical capture cross section ath is given by 
the product of the DC cross section ctdc and the spec- 
troscopic factor C^S of the final state: 



fth = C S X (Juc 



(2.1) 
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The DC cross section depends on the square of the over- 
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lap between the scattering wave function XLi,Jiir), the 
bound state wave function ujqj^^ j^ir)^ and the electro- 
magnetic operator qEi.E2,mi j^^^ j^2^ and Ml tran- 
sitions: 



fDC 



E1.E2,AI1 



UN.Lj,jj{r)dr (2.2) 



where L^, Ji and Lf,Jf are the angular momenta in the 
initial scattering wave function xi''') ^-nd the final bound 
state wave function u{r). N is the number of nodes in the 
bound state wave function that takes into account the 
Wildermuth condition. The total number of oscillator 
quanta Q — 2N + Lf = 3 for three nucleons in the Ip 
shell leads to two cluster states in ^Be. There is a first 
state with one node {N = 1) and angular momentum 
L = 1, and a second state without node {N = 0) and 
L = 3. Both states are split into dubletts because of the 
spin 5 = 1/2 of the "^He nucleus. Both L = 1 states 
are located below the "^He-a threshold in ^Be, whereas 
the L = 3 states are located above the threshold and 
may appear as resonances in the ^He(a,7)^Be capture 
reaction. Properties of the cluster states are given in 
Table HI and a simplified level scheme of ^Be is shown in 
Fig. [2 



TABLE I: Level scheme of '^Be. Experimental data are taken 
from [13]. The potential strength parameters are adjusted to 
reproduce the energies of the bound and quasi-bound states 
(see text). 



L J" 



E* 

(keV) 



cxp — -S'calc 

(keV) 



-L cxp 

(keV) 



Tcalc 

(keV) 



1 3/2" 
1 1/2" 429.1±0.1 
3 7/2' 4570 ± 50 
3 5/2" 6730 ±100 



-1586.6" (Ti/a = 53.22 d) 
-1157.5" (r = 192 ± 25 fs) 
2983.4" 175 ± 7 
5143.4" « 1200 



145 
i 1100 



"potential adjusted 

The bound state wave functions u{r) are shown for 
the L = 1, J"" — 3/2" ground state and the L — 1, 
= 1/2- first excited state in '^Be in Fig. H Fig. [3] 
shows the integrand of Eq. (|2.2p . i.e. the overlap of the 
scattering wave function x(^); the electromagnetic oper- 
ator O, and the bound state wave function u{r) for dif- 
ferent transitions in the •^He(a,7)^Be capture reaction at 
a very low (100 keV) and at a higher energy (3MeV). A 
detailed discussion will be given in the following Sec. IIIII 

The basic ingredient for the calculation of the DC in- 
tegral in Eq. (|2.2p are the potentials for the entrance 
channel (elastic scattering) and the exit channel (bound 
state wave function). As soon as the potentials are fixed, 
the DC integral is calculated without further adjust- 
ment of parameters to experimental capture data. The 
spectroscopic factor C^S of the bound state is taken - 
exactly as in |31j - from theory in the present study: 
C'^S{3/2-) = 1.174 and C^S{l/2-) = 1.175 g^. The 
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FIG. 1: Simplified level scheme of ^Be with the L — 1 and 
L = 3 cluster states. Excitation energies E* are given in keV. 
All data are taken from 14711. 



spectroscopic factor is considered as an absolute normal- 
ization for the calculated capture cross section, similar 
to the procedure in [s^]- It is beyond the scope of the 
present paper to discuss the relation between the asymp- 
totic normalization coefficient and the spectroscopic fac- 
tor as e.g. in [5l| because the main conclusions of this 
work are not affected. 

The potential for the "^He — a system is the sum of the 
central nuclear potential, the spin-orbit potential, and 
the Coulomb potential. The real part of the central nu- 
clear potential is calculated from the double folding pro- 
cedure [H, nil that is scaled by a strength parameter A 
which is about 1.4 — 1.8 in this study. The underlying 
nuclear densities are derived from the measured charge 
density distributions ^54i | . The imaginary part of the nu- 
clear potential can be set to zero because there are no 
open channels below 4MeV except the relatively weak 
■^116(0!, 7)^Be capture. The spin-orbit potential is taken 
in the usual Thomas form proportional to 1/r x dV/dr, 
again scaled by a spin-orbit strength parameter Xls- Fi- 
nally, the Coulomb potential Vc is calculated from the 
homogeneous charged sphere model with a Coulomb ra- 
dius Rc identical to the root-mean-square radius of the 
folding potential: Rc — r^^s — 2.991 fm. Further details 
can be found in [31]. The total potential is given by: 



V{r) = A Vpir) + Xls 

r ar 



LS + Vcir) (2.3) 



with the unsealed (A = 1) folding potential Vpir). The 
factor fm^ in the spin-orbit potential is added to obtain 
a dimensionless strength parameter Xls- L, S, and J — 
L + S are the orbital, spin, and total angular momenta 
in units of h. 
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FIG. 2: Bound state wave functions u(r) for the L = 1, 
J'^ = 3/2^ ground state and the L ^ 1, T = 1/2" first 
excited state in '^Be in linear scale (upper) and v?(r) in loga- 
rithmic scale (lower). The wave functions are very similar in 
the nuclear interior. The slope in the exterior reflects the dif- 
ferent binding energies [Eb = — 1587keV for the 3/2~ ground 
state and Eb = — 1158keV for the l/2~ first excited state). 



One main advantage of folding potentials is the small 
number of adjustable parameters. The shape of the 
potential is fixed by the folding procedure. Only the 
strength parameters A and Als have to be adjusted to 
experimental data. Obviously, the small number of ad- 
justable parameters improves the predictive power of the 
calculations. 

The shape of the folding potential for ■^He — a is almost 
Gaussian. Consequently, the spin-orbit part has almost 
the same shape, and the sum of central and spin-orbit 
potential again has the same shape. In other words, the 
spin-orbit potential slightly increases the potential for 
J = L -1-1/2 waves and slightly decreases the potential for 
the J = L~\/2 waves but keeps the shape of the folding 
potential. Alternatively, the same effect can be obtained 
if the spin-orbit potential is set to zero and instead the 
strength parameter of the central potential becomes J- 
dependent. This approach has been followed e. g. in [s^ : 
there very similar results to our previous work 31| were 
obtained using empirical Gaussian potentials with a J- 
dependent depth Vq. 



III. GENERAL REMARKS 

The nucleus ^Be and the '^IIe(Q;,Q;)"^He elastic scatter- 
ing and '^IIe(Q;,7)^Be capture reactions are textbook ex- 
amples for the successful application of a simple two-body 
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FIG. 3; Overlap of the scattering wave function x('')i electro- 
magnetic transition operator O, and bound state wave func- 
tion ti(r) for different transitions of the '^He(a,7)^Be cap- 
ture reaction. Upper: Non-resonant E\ s-wave capture to 
the 3/2~ ground state and l/2~ first excited state at ener- 
gies E = 0.1 MeV and 3.0 MeV. Middle: Non-resonant El 
d-wave capture to the 3/2~ ground state at iJ = 0.1 MeV and 
3.0 MeV. Lower: Non-resonant and resonant E2 /-wave cap- 
ture to the 3/2" ground state aX E = 0.1 MeV and 3.0 MeV. 
Further discussion see text. 



model because of the strong internal binding energies of 
the two constituents ^He and ''He. Early calculations 
have been performed for the '^He(a,7)^Be reaction using 
hard-sphere phase shifts for the entrance channel, thus 
considering external capture only. Already these early 
calculations have successfully reproduced the cross sec- 
tion of the '^He(Q;,7)^Be capture reaction at low energies 
[20I . [21I [2^ . Nowadays it must be the aim of theoreti- 
cal studies to reproduce simultaneously the •^He(Q;,Q;)'^He 
elastic scattering cross sections and phase shifts, the 
■^He(Q!,7)^Be capture reaction, and electromagnetic prop- 
erties of the ""Be nucleus. Often also the mirror nucleus 
'^Li and the mirror reaction ^H(Q!,7)^Li are considered. 
The present study focuses on the new results for the 
■^He(a,7)^Be reaction [1]; the ^Li mirror system was al- 
ready studied in our earlier work jsij . 
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Only very few radiation widths or lifetimes r have 
been measured for the ^Be nucleus. The first excited 
l/2~ state in ^Be decays by a Ml transition to the ^Be 



192 ± 25 fs 



ground state with a lifetime of = 
responding to B(Ml,3/2- ^ U^' 
Following the formalism in [ssl . Issl l56l | , the reduced tran- 
sition strength B{M1, 3/2~ l/2~) is given by 



47[ cor- 
1.87 ± 0.25 Aii,. 



B(M1) 



An 



\< "3/2-|wi/2 



- > 



(3.1) 

with the effective cluster-cluster orbital gyromagnetic 

factor G = {ZiAj + ZjA?)/ [^1^2(^1 + A2)] = 0.5952 

for ^Bc = ^Re (g) a and ^(f He) = -2.1276 As 

already pointed out in [30|, the overlap of the bound 

I 1 2 

state wave functions |< M3/2-IM1/2- >| is close to 

unity because both wave functions u{r) are very simi- 
lar. This similarity is also obvious for the wave functions 
in this study, see Fig. [21 here the overlap deviates by 
less than one per cent from unity. Consequently, the re- 
duced transition strength B{M1) is practically defined 

by the factor ^ [2fi{^Re) ~ fiNG]^ = 1.873 /z^ which 
is in excellent agreement with the experimental result of 
1.87±0.25^^. The excellent reproduction of the B{M1) 
transition strength between the 3/2^ ground state and 
the l/2~ first excited state in ''Be is thus a more or less 
trivial result for any two-body calculation with realistic 
bound state wave functions. 

The E2 contribution of this transition is orders of mag- 
nitude smaller than the Ml transition because of the 
relatively low transition energy of — 429 keV and 
the strong dependence of E2 transitions: one E2 
Weisskopf unit corresponds to a radiation width of about 
r^(i?2) « lOneV which has to be compared to the ex- 
perimental Ml width of = 3.4 meV. 

No further experimental data for radiation widths in 
^Be exist in 47 1. However, the recent experiment has 
measured the cross section of the ^He(a,7)'^Be reaction 
in the 7/2~ resonance for the first time, and it was pos- 
sible to derive the radiation width F^.o in [2]. It will 
be shown in Sec. lIVBl that the DC model is able to re- 
produce the experimental value within its relatively large 
experimental uncertainty. Contrary to the above studied 
Ml transition, this result is not trivial because it requires 
the overlap of the 7/2~ scattering wave function, the E2 
operator, and the 3/2~ ground state wave function (see 
Fig. [21 lower part). 

Often, "direct capture" is also named "extranuclear 
capture" or "non-resonant capture" . The origin of these 
somewhat misleading names comes from the early cal- 
culations by Christy and Duck [l^l and Tombrello and 
coworkers [2l|, A good description of the experi- 

mental ^He(a,7)^Be capture cross section was obtained 
in [21I, I23] using hard-sphere phase shifts for the 
non-resonant s- and d-waves (that are by definition non- 
resonant). As can be seen from Fig. [31 (upper and mid- 
dle parts), the main contribution of the El transitions 
from incoming s-waves and d-waves to the bound L = 1 



states comes indeed from the nuclear exterior, thus val- 
idating the hard-sphere approximation for the s-waves 
and d- waves in (20II21I [2^ . However, as can be seen from 
Fig. [31 (lower part), the main contribution of the E2 tran- 
sition from the incoming /-wave to the 3/2^ ground state 
comes from the exterior at the low energy oi E ^ 100 keV, 
but is shifted to the interior at E = 3.0 MeV which is in 
the 7/2~ resonance. Thus, the terminus "direct capture" 
is neither identical to "extranuclear capture" nor to "non- 
resonant capture" . This has also been illustrated for the 
3+ resonance in the ^H(a,7)^Li reaction in [5^. 



IV. RESULTS 
A. Elastic scattering cross sections and phase shifts 

As already pointed out above, the parameters of the 
potential - i.e. the potential strength parameters A and 
Als - have to be adjusted to experimental data. Elastic 
scattering phase shifts are typically used for this adjust- 
ment. Such phase shifts have been derived from angular 
distributions and excitation functions over a broad range 
of energies HH, SI, SI, see Fig. [2 

The data of Spiger and Tombrello cover the energy 
range slightly below the 7/2~ resonance and range from 
about 2.5 MeV up to 10 MeV [4l|. The data of Boykin 
et al. cover the energy range from about 1.9 MeV up 
to 4 MeV. However, no numerical results are given for 
the important d-waves. Instead, it is pointed out that 
"the values obtained at the various energies were scat- 
tered with no discernible trend in a band between —4° 
and -|-4°" [i^l. The data from Hardy et al. start above 
the 7/2~ resonance and range from about 3.3 MeV up to 
7.7 MeV Q. 

The adjustment of the potential strength parameters 
A and Als is done as follows. The procedure is almost 
identical to the previous work [Slj with a modification 
of the spin-orbit potential. First, the strength param- 
eter A is adjusted to the s-wave phase shift. This can 
be done unambiguously because the s-wave is not af- 
fected by the spin-orbit potential. An excellent agree- 
ment with the experimental phase shifts is obtained with 
A(L = 0) = 1.452 (full black line in Fig. [H upper part). 
Additionally, the results are shown for a variation of the 
strength parameter A between 1.4 and 1.6 (green dashed, 
red dotted, and blue narrow-dashed lines) that will be 
important for the analysis of the '^He(Q!,7)^Be capture 
cross section in Sec. IIVBI 

The parameter X{L = 0) = 1.452 is also used for the 
d-wave phase shifts. As can be seen from Fig. [4l there is 
excellent agreement with the experimental d-wave phase 
shifts over a broad energy range. There is no evidence for 
a different behavior of the d3/2 and d5/2 phase shifts from 
the data in [4l], |4|. Thus, the spin-orbit potential 
has to vanish for the d- waves, and \ls{L — 2) — 0. It 
is interesting to note that any deviation of the potential 
strength from A = 1.452 leads to deviations of the calcu- 
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FIG. 4: (Color online) Phase shifts for elastic ^He(a,Q)^He 
scattering. Experimental data are taken from [Zll . |43 . |4^ . 
Calculations are performed using the potential strength pa- 
rameters A''^™ = 1.452, Ar™ = 0, A""*"* = 1.830 and 
>^°LS = -0.173 (standard case, full black line). The parame- 
ter A^^"' is varied between 1.40 and 1.60 (A°™" = 1.40: green 
dashed line; A'^^™ = 1.50: red dotted line; A''^"" = 1.60: blue 
narrow-dashed line), and the spin-orbit strength of the odd 
partial waves has also been used for the even partial waves 
(A''™" = 1.452, Af^^l" = -0.173: magenta dash-dotted line). 
Further discussion see text. 



lated d-wave phase shifts from the experimentally vanish- 
ing data. Combining the results for the s-wave and the 
d-waves, the potential strength parameters for even par- 
tial waves are A''™" = 1.452 and A^^'J" = 0. A test with 
the spin-orbit strength derived from the /-waves (see be- 
low) also shows clear disagreement to the experimental 
d-wave phase shifts (magenta dash-dotted line in Fig. [5]). 

A noticeably higher value for the potential strength A 
is found for the odd partial waves. A"^'* and have 
been adjusted to reproduce the energies of the /-wave 
resonances (see Fig. S]). This leads to A°'^'^ = 1.830 and 
— —0.173. Using these values, an excellent agree- 
ment for the /-wave phase shifts is obtained in both non- 
resonant and resonant energy regions. In particular, the 
resonance widths are properly reproduced (see also Table 
m . Good agreement with the experimental p-wave phase 
shifts is also obtained with the above A°'^^ and X°^g. The 
parameters A"'*'^ and X°^g were not varied in the follow- 
ing study because any change of these parameters shifts 
the /-wave resonances dramatically in energy. 

The new experimental ■^IIe(a,7)^Be capture data cover 
an energy range that is at the lower end of the experi- 
mental phase shift data of [l!], [42, El] . Fortunately, fur- 
ther experimental angular distributions of ^He(a,a)'^He 
elastic scattering are available in [3l|, IH, in the rele- 
vant energy range of the new ■^IIe(Q;,7)^Be experiment 0]. 
These experimental angular distributions are compared 
to the calculated cross sections using the above deter- 
mined parameters A'=™" = 1.452, Af|" = 0, A"'^'^ = 1.830 
and Xf^ = -0.173. The result is shown in Fig. [5] (fuU 
black lines). Additionally, again the sensitivity on the po- 
tential strength parameters A°'^°" and A|^™ is analyzed 
(colored dashed and dotted lines in Fig. [5]). 

The conclusions from Fig. [5] are clear. As expected, 
an excellent agreement between the experimental data 
and the calculations is only obtained using the same pa- 
rameters as derived from the phase shifts of [4l|, 
in Fig. [31 From the shown angular distributions the pa- 
rameter A''™" may vary between 1.40 and 1.45; a care- 
ful inspection of the phase shift calculations in Fig. [4] 
shows that the low-energy s-wave phase shifts of [i^ and 
the overall energy dependence are best described using 
^even _ 1.45; howcvcr, the data from [4lj are better 
desribed using A°™" — 1.40 (green dashed hue in Fig. [4]). 
Because of the better overall description of the s-wave 
phase shifts and the excellent description of the d-wave 
phase shifts I take A'=™" = 1.45, Af|" = as the stan- 
dard calculation in the following, keeping A"^*^ = 1.830 
and A°1j^ = -0.173 fixed for the odd partial waves. An 
additional calculation with A'^™" = 1.45, A^^^" = -0.173 
(instead of X f^ = 0) is almost identical to the previous 
calculation in [3l| except minor technical modifications 
(e.g. larger integration range and smaller step size). 

A minor disagreement between the calculated angular 
distribution and the experimental data can be seen at 
E — 2715 keV. This is not very surprising because this en- 
ergy is at the low-energy tail of the 7/2~ resonance. The 
width of this resonance is slightly underestimated: the 
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FIG. 5: (Color online) Angular distributions of ^He(Q,Q)^He 
elastic scattering, normalized to Rutherford scattering. Ex- 
perimental data are taken from [3ll . |43 ]. The linestyle and 
color coding is identical to Fig. |4] the dotted black line at 
E — 2715 keV shows the artificial enhancement of the contri- 
bution of the 7/2~ resonance (see text). 



calculated width is Tr 



145 keV which has to be com- 



pared to the experimental width of 



175±7keV 



(see also Table IJ). A much better agreement is obtained 
if the resonance contribution is artificially enhanced in 
the calculation by an increase of the phase shift of the 
7/2~ partial wave by about 5° from 10° in the standard 
calculation to 15° (dotted black hne in Fig. [5]). The de- 
viation is clearly not related to the potential strength of 
the even partial waves which will be most important for 
the analysis of the •^He(Q!,7)^Be capture cross section in 



the next Sec. HVB] 

Summarizing the above, an excellent agreement be- 
tween all experimental scattering data (phase shifts and 
angular distributions) and the calculated values is ob- 
tained using the above derived potential strength pa- 
rameters A'=™" = 1.452, A^;'^" = 0, A°'i'i = 1.830 and 
A^lj* = —0.173 (standard case). The results of all scat- 
tering experiments are consistent with each other and can 
be well described using the present model. 

The scattering wave functions can now be calculated 
without further adjustment of parameters. Because of 
the dominating external capture (see Fig. [3]), good agree- 
ment between the predicted and experimental capture 
cross sections has to be expected. This will be further 
analyzed in the next Sec. IIVBI 



B. The ^He(Q,7)^Be capture cross section 

In addition to the scattering wave function the 
determination of the DC cross section in Eqs. (12. ip and 
(|2.2|) requires the calculation of the bound state wave 
function u{r). Because of the contributions from the 
nuclear exterior (see Fig. [3]), it is essential that the 
calculated bound state wave function has the correct 
asymptotic behavior that is defined by the binding en- 
ergy. This leads to a slight readjustment of the potential 
strength parameter A for the bound states compared to 
the above results for the scattering phase shifts. The re- 
sults are A(3/2~) = 1.836 for the ground state of ''Be 
at = -1587keV and A(l/2") = 1.799 for the first ex- 
cited state of ^Be at E ^ -1158 keV {E* = 429 keV). No 
spin-orbit potential has been used here. 

The remaining calculation of the DC cross section is 
straightforward and does not require any adjustment to 
experimental capture data. The results are shown in 
Fig. [6] with a comparison to the available experimental 
data and in Fig. [7] in a broader energy range. The cal- 
culations take into account all possible El, E2, and Ml 
contributions from the incoming scattering waves with 
L = — 3. However, as has been shown in several pre- 
vious studies, the El contribution is dominating at all 
energies. Only in the 7/2" and 5/2~ resonances a no- 
ticeable E2 contribution is found, and the Ml contribu- 
tion remains negligible over the whole energy range. A 
decomposition into the various contributions has alr eady 
been shown in Fig. 5 of i2a]. Fig. 8 of |2a], Fig. 8 of [2|, 
Fig. 2 of [Hi, Fig- 7 of [3lf, Fig. 2 of Fig. 8 of 
and Fig. 5 of ^8], and is thus not repeated here. 



The standard calculation (using the potential strength 
parameters from the adjustment to the experimental 
phase shifts) shows an excellent agreement with the total 
S-factor and the branching ratio R — cr429/cro in the low- 
energy region below 1 MeV. At higher energies the calcu- 
lation agrees with the old data of Parker and Kavanagh 
[3| , but is significantly lower than the recent experimental 
data of 

The sensitivity of the ■^He(Q;,7)^Be capture cross sec- 
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detailed discussion will be given in Sec. [V] 
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FIG. 6: (Color online) Astrophysical S-factor and branching 
ratio R — cr429/cro for the ^He(a,7)^Be capture reaction. Ex- 
perimental data are taken from 0, 0, i, @, Q, H i, [13, EH, 
The llnestyles and color codes are the 
same as in Figs. |3] and [S] The vertical arrows on top indi- 
cate the energies of the angular distributions in the previous 
Fig. [5] where the elastic scattering cross sections are well re- 
produced. These scattering data cover the energy range of the 
recent ^He(a,7)^Be capture data 2] (shown as open squares). 



tion on the underlying potential has been studied in the 
same way as in the previous Sec. IIV A! bv a variation of 
the A and X^s potential strength parameters. The results 
are also shown in Figs. [S] and [71 In general, an increas- 
ing potential strength leads to an increasing S-factor. A 




4000 
E (keV) 



8000 



FIG. 7: (Color online) Same as Fig. (6] but with an extended 
energy range up to 8 MeV. For better visibility of the calcu- 
lated curves, the experimental data have been omitted. The 
linestyles and color codes are the same as in the previous 
Figs.HISl andU 

Before entering the detailed discussion of the shown 
■^He(Q!,7) ''Be capture cross sections, the presented results 
are completed by an analysis of the two L = 3 resonances 
and a study of the low-energy behavior of the astrophys- 
ical S-factor. For both resonances it is obvious that the 
total width F is practically identical to the alpha particle 
width Fq,. As already shown in Table [H the calculated 
widths Fcaic agree reasonably well with the experimental 
results Fexp- This is also obvious from correct description 
of the resonant behavior of the L = 3 elastic phase shifts 
in Fig. 12 

The 7/2~ resonance at 2983 keV decays by an E2 tran- 
sition to the 3/2^ ground state of ^Be whereas no E2 
transition to the first ecited l/2~ state is possible. A 
Breit-Wigncr fit to the calculated capture cross section 
in this resonance leads to a radiation width F-, = F-^.o = 
37meV and a resonance strength oj^ = 150 meV. This 
prediction is in reasonable agreement with the experi- 
mental result of Lu^ = 330 ± 210 meV [2]. 

The 5/2~ resonance at 5143 keV was not measured 
yet in the ^He(Q;,7)''Be reaction. The DC model pre- 



dicts F 



7,0 



50meV for the E2 transition to the 3/2 



ground state and F 



7,429 



185 meV for the E2 tran- 
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sition 

r. 



the first excited 1/2 state. This leads to 
235 meV or a total resonance 



- 7,429 



7 ^ 7,0 

strength ui-f w 0.7 eV. Because of the huge width of this 
resonance, any determination of resonance parameters 
has relatively large uncertainties. An additional theo- 
retical uncertainty is related to the opening of the ^Li-p 
channel slightly above 4 MeV. 

It is interesting to note that the different decay pat- 
terns of the 7/2~ and 5/2^ resonances are clearly visible 
in the branching ratio R although the E2 contribution to 
the total cross section is small compared to the dominat- 
ing El transitions. Fig. [7] clearly shows a smaller branch- 
ing ratio in the 7/2~ resonance around 3 MeV which can 
only decay to the 3/2^ ground state, and a larger branch- 
ing ratio in the 5/2~ resonance around 5 MeV which pre- 
dominantly decays to the l/2~ first excited state. 

The S-factor at very low energies is extrapolated down 
to 5(0) using the following procedure. The calculated 
cross sections are fitted using a second-order polynomial 
up to 500 keV. The results for S{0) and 5'(0)/S'(0) are 
listed in Table [IT] for the different potentials shown in 
Figs. H El El and in The standard case leads to 5(0) = 
0.53 keVb and 5'(0)/5(0) = -0.73 MeV'^, in agreement 
with the compilation of solar fusion cross sections by 
Adelberger et al. 18,] who recommend 5(0) = 0.53 ± 
0.05keVb and 5'(0)/5(0) = -0.57MeV-\ and the 
NACRE compilation [13 where 5(0) = 0.54 ± 0.09 keVb 
and 5'(0)/5(0) = -0.96 MeV"! are gjven. The re- 
cent summary of Cyburt and Davids [40| recommends 
a slightly higher value of 5(0) = 0.580 ± 0.043 keVb 
and 5'(0)/5(0) = -0.92 ± 0.18MeV-i. Variations of 
the order of about 10% for the slope of the S-factor 
at zero energy are found using the different potentials. 
Such variations are also Wmcal for the different theoret- 
ical ^lculation_ynj2^j21^^ [H, IH, 




34 |35|, [361, IM 1391 . The adopted 
40| show somewhat larger variations. 



|29|, |30|, IM 13 

5'(0)/5(0) [1 

The average of the three compilations of 5'(0)/5(0) 
— 0.82MeV~^ is in reasonable agreement with the calcu- 
lated 5'(0)/5(0) = -0.73MeV^i in the standard case 
and in perfect agreement with the calculation using 
js^cvcn ^ 2^ 4Q 5"(0)/5'(0) = -0.79MeV-i is found. 

However, larger values A°™" > 1.50 correspond to lower 
5'(0)/5(0). 

Using the standard case potential, 5(0) = 0.53 keVb 
is predicted from the present study. This prediction de- 
pends on the theoretical spectroscopic factors C^5 of the 
bound state wave functions. From the excellent agree- 
ment between the experimental and calculated capture 
cross sections and branching ratios at low energies (see 
Fig. E]) it is obvious that the chosen theoretical spectro- 
scopic factors 2^ are correct within about 10%. Any 
minor readjustment of these spectroscopic factors will af- 
fect the the calculated capture cross sections linearly at 
all energies, but will not affect the calculated energy de- 
pendence of the '^He(Q;,7)''Be capture cross section. Also 
the normalized slope 5'(0)/5(0) will not be affected by 
a readjustment of the spectroscopic factors. 



TABLE II; Extrapolated astrophysical S-factor S{Q) at zero 
energy and its normalized derivative 5''(0)/S'(0) for diflFerent 
potentials. The strength parameters A°'*'* = 1.830 and = 
—0.173 have not been varied. For comparison, the results of 
recent compilations [isl . [Tgl . |40| | are also listed. 



^cvcn 


\ even 


S(0) 
(keVb) 


S'(0)/S'(0) remarks 
(MeV-^) 


1.45 


0.0 


0.530 


-0.731 


standard case 


1.40 


0.0 


0.497 


-0.791 




1.50 


0.0 


0.557 


-0.666 




1.60 


0.0 


0.611 


-0.536 




1.45 


-0.173 


0.530 


-0.733 


close to Ref. [3^1] 






0.53 ±0.05 
0.54 ±0.09 


-0.57 
-0.96 


Adelberger et al. [18] 
NACRE [19] 






0.580 ± 0.043 


-0.92 


Cyburt and Davids [40] 



The chosen DC model does not require effective 
charges and uses q — +2e for the involved ^He and ^He 
nuclei. An additional introduction of effective charges in 
the model may affect the absolute values of the DC cross 
sections. Because of the dominating El contribution to 
the capture cross section, the introduction of effective 
charges does not affect the calculated energy dependence 
in a significant way. Similar to the statements on the 
spectroscopic factor in the previous paragraph, effective 
charges do also not affect the normalized slope 5'(0) /5(0) 
of the astrophysical S-factor. Although effective charges 
cannot be determined from scattering phase shifts, it is 
obvious from the excellent agreement between the ex- 
perimental and calculated cross sections at low energies 
that the effective charges cannot deviate strongly from 
the used real charges. 



DISCUSSION AND CONCLUSIONS 



It has been shown in the previous Sects. [TVH and 
lIVBI that the standard calculation with the parameters 

Xc^cn ^ ^cvcn ^ ^odd ^ gg^ ^odd ^ _0.173 

is able to describe elastic scattering phase shifts for 
'^IIe(a,a)'^He reaction in the full energy range up to about 
8 MeV and the '^He(a,7)^Be capture cross sections to 
the ^Be ground state and first excited state at low en- 
ergies. At higher energies the calculation agrees with 
the old data [4| but disagrees with the new data 0. In 
the following discussion the theoretical uncertainties will 
be analyzed, and the question will be answered whether 
the new data of can be described within the present 
model. Most of the considerations will also be valid for 
other models because of the dominating external contri- 
butions to the ■^IIe(a,7)^Be capture cross section. The 
discussion will focus on the dominating El contributions 
from the incoming s-wave and d-wave. 

The analysis of elastic scattering clearly shows that the 
s-wave and d-wave phase shifts can only be reproduced 
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using the potential strength of the standard calculation 
(see Figs. Hand [5]) . For the s-wave A°™" = 1.40 - 1.45 
is the acceptable range. The d-wave phase shifts which 
remain close to zero up to about 8MeV can only be de- 
scribed with A values very close to A"™" = 1.45. A value 
of A°™" > 1.50 can clearly be excluded. This is further 
strengthened by the analysis of the slope S''(0)/S'(0) at 
low energies, see Table HIl 

The integrand of the DC integral in Eq. (|2.2p . i.e. the 
overlap between scattering wave function electro- 
magnetic operator O^^, and bound state wave function 
u{r), is shown for the standard calculation in Fig. [3] The 
upper part shows the El transition from the s-wave to 
the ground state and first excited state, and the mid- 
dle part shows the El transition from the d-wave to the 
ground state in ^Be. As usual, the integrand is shifted 
towards the exterior for lower energies. But it has to 
be pointed out that all curves are clearly dominated by 
contributions from far outside the nucleus at both con- 
sidered energies E = 100 keV and 3MeV; i.e., also at 
the high energy E = 3MeV the external contribution 
dominates. Thus, the result of the DC calculation is es- 
sentially defined by the behavior of the scattering wave 
function x(r) and the bound state wave function u{r) at 
radii between 5 fm and 35 fm that is much larger than the 
size of ^Be. Any calculation with realistic nuclear poten- 
tials should provide a more or less similar capture cross 
section if the potential is able to reproduce the elastic 
phase shifts. This is particularly true for the calculated 
energy dependence. 

The absolute value of the bound state wave function 
u{r) in the external region depends indirectly on the cho- 
sen potential shape. The slope of the external part of 
the bound state wave function u(r) is well defined by the 
binding energy, but the normalization of the bound state 
wave function u(r) 

/>oo 

/ u^{r)dr = l (5.1) 

JO 

has its main contribution from smaller radii r < 5 fm, 
i.e. from the nuclear interior and surface (see Fig. [5]). In 
line with the considerations in [HO], the chosen spectro- 
scopic factors C'^S « 1.2 for the ground state and first 
excited state in ^Be simply mean here that the bound 
state wave functions u{r) in the nuclear exterior have to 
be scaled by a factor yCPS « 1.1 to bring the calculated 
DC cross section into agreement with the experimental 
data at low energies. The choice of a different poten- 
tial leads to a different shape of the wave function u(r) 
and its asymptotic behavior that can be compensated by 
a different spectroscopic factor. Finally, the theoretical 
cross section, in particular its energy dependence, will re- 
main almost unaffected - provided that the phase shifts 
are well reproduced. 

Nevertheless, an attempt has been made to repro- 
duce the new '^He(a,7)^Be capture data within the 
present model. In general, an increasing nuclear poten- 
tial strength decreases the Coulomb barrier and thus in- 



creases the calculated capture cross section. A value of 
^ovcn ^ is required to reproduce the larger cross sec- 
tion data of [2]. The calculated capture cross section 
at low energies is increased by about 15% compared to 
the standard case which is slightly higher than the al- 
lowed range of the precision experiments at low energies. 
This larger potential strength also leads to a strong dis- 
agreement with the scattering data. In particular, the 
elastic scattering angular distributions in Fig. [5] cannot 
be described with the increased value of A"™" = 1.6. 
These angular distributions have been measured in an 
energy range that covers the energy range of the recent 
■^He(Q!,7)'^Be capture data (indicated by the arrows 
on top of Fig. [6]). Taking into account the above consid- 
erations on the radial behavior of the overlap integral in 
Eq. (|2.2p . it can be excluded that the DC model is able 
to reproduce simultaneously ■^He(a,a)'^He elastic scatter- 
ing SH, m S, m and 3He(a,7)^Be capture at 
energies around 2 — 3 MeV. 

It is interesting to note that the branching ratio R = 
CT429 /co is practically not affected by the variation of the 
potential strength of the s- and d-waves because both 
transitions are increased in the same way when A''™" is 
increased from its standard value of 1.45 to about 1.6 to 
fit the higher capture data of [2^]. However, the branching 
ratio R is significantly affected by the spin-orbit poten- 
tial of the d-wave. Using the spin-orbit potential strength 
from the odd partial waves also for the even partial waves 
(A^;^™ = -0.173 instead of the standard case A^^^" = 0), 
the effective strength of the potential for the ^5/2 (^^3/2) 
wave is increased (decreased). As pointed out above, an 
increased potential strength leads to increased capture 
cross sections. Because there is only an El transition 
from the wave to the 3/2^ ground state, but no 
El transition to the l/2~ first excited state, the ground 
state contribution is enhanced, and thus the branching 
ratio R decreases at higher energies (see Figs. in]and[71). 
However, at low energies the dominating contribution to 
the capture cross section comes from the s-wave which is 
not affected by the spin-orbit potential, and the branch- 
ing ratio R at very low energies is practically identical to 
the standard case. 

It is difficult to find an explanation for the discrep- 
ancy between the calculated '^He(a,7)''Be capture cross 
sections and the recent experimental data [4]. Experi- 
mental scattering data have been obtained in a series of 
independent experiments [Ml, |4l|, IH, [Ulil, SI El] that 
agree well with each other. The determination of the 
scattering wave function in the DC integral in Eq. (j2.2|) 
is thus very well-defined. Although the requested accu- 
racy of a few per cent for astrophysical modelling has 
not yet been achieved for the '^He(a,7)^Be capture reac- 
tion, all experimental capture data at low energies agree 
within each other at a 10 — 15 % level. The spectroscopic 
factors - taken as normalization factors in Eq. (|2.2p - 
are thus also well-defined at this level. The same level of 
accuracy should then be expected for the energy range 
of the new '^He(a,7)^Be capture data up to about 
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3 MeV. Instead, the standard calculation underestimates 
the new capture data by about 30 — 40%. It has to be 
pointed out that a wrong normalization of the new data 
[2] is also very unlikely because the experiment has stud- 
ied the "^116(0, 7)^Be capture reaction very carefully using 
three different independent techniques (recoil detection, 
7 detection, and activation) that agree very well with 
each other at the level of a few per cent. Of course, an 
independent confirmation of the new data in is nev- 
ertheless highly desireable, in particular, because of the 
discrepancy between the old data in Q and the new data 
in 0. 

The only remaining explanation is that the DC model 
itself is not appropriate for the description of the 
■^IIe(ck,7)^Be capture reaction. However, even "exotic" 
theoretical solutions like other electromagnetic transi- 
tions (except El, E2, and A/1) or further unobserved 
bound states in ^Be can be excluded at the requested 
level of 30 — 40 % from the agreement of the different 
experimental techniques in Q . And because of the dom- 
inating exteral contributions to the DC cross section, 
any theoretical calculation that reproduces the scattering 
phase shifts, should provide a similar theoretical energy 
dependence and thus underestimate the new capture data 
for the '^He(a,7)''Be reaction in 0] in a similar way. 



model together with folding potentials. It is shown that 
the elastic scattering phase shifts are well reproduced 
over a broad energy range. Whereas the low-energy 
capture data are well described simultaneously with the 
phase shifts, there is a significant underestimation of the 
recent capture data , but agreement with the old data 
of [3| . As the DC calculation is very well constrained by 
experimental scattering data over the whole energy range 
studied in , it is difficult to find an explanation for the 
discrepancy between the calculated and the recently mea- 
sured capture cross section of the ■^He(a,7)^Be reaction. 
Further eperimental and theoretical effort is needed to 
resolve this surprising problem. 
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